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MOORE-GIBSON-THOMPSON EQUATION WITH 
MEMORY, PART II: GENERAL DECAY OF ENERGY 

IRENA LASIECKA AND XIAOJUN WANG* 


Abstract. We study a temporally third order (Moore-Gibson-Thompson) 
equation with a memory term. Previously it is known that, in non- 
critical regime, the global solutions exist and the energy functionals de¬ 
cay to zero. More precisely, it is known that the energy has exponential 
decay if the memory kernel decays exponentially. The current work is a 
generalization of the previous one (Part I) in that it allows the memory 
kernel to be more general and shows that the energy decays the same 
way as the memory kernel does, exponentially or not. 


1. Introduction 


We study the energy decay of Moore-Gibson-Thompson(MGT) equation 
with a viscoelastic term 


( 1 ) TUttt + (xutt + c^Au + bAut 

with initial data 



s)Au{s)ds = 0, 


(2) u(0) = no, tit(O) = ni, ntt(O) = U2, 

where r, c, b are parameters inherited from modeling process, see m and 
references therein. The constant a can be scaled out; we keep it, however, 
for notational consistency with [12]. A is a positive self-adjoint operator 
defined in a real Hilbert space R. The convolution term g{t — s)Au{s)ds 
reflects the memory effect of viscoelastic materials; the “memory kernel” 
g{t) : [0, oo) —>• [0, oo) directly relates to whether or how the energy decays. 
Without this memory term, it is known the MGT equation has exponential 
energy decay in the non-critical regime, where 7 = 0 — ^ > 0, see m- 
In an earlier work |18] . we studied ([T|) with a nontrivial g{t), but focusing 
on the case where g{t) has exponential decay. We were able to get expo¬ 
nential decay of the energy for three types of memories in the non-critical 
regime. Here we study the case where the memory kernel has a more general 
decay rate. For sake of clarity, in this work we restrict our attention to one 
of the three types of memories introduced in |18j . 

Notations: 

• g{t): memory kernel. 

• G{t) = fg g{s)ds: strength of memory. 

• H: real Hilbert space. 
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• 11 • 11: norms on H. 

• g o h = g{t - s)||/i(t) - h{s)\\‘^ds, g{-) € C{R^),h{-) G H. 

1.1. Main results. Our work shows that the energy decay rate of system 
where memory effects get involved, is determined solely by the memory 

kernel: if g{t) decays exponentially, then the energy decays exponentially 
too; if g{t) decays slower, then the energy decays slower as well. 

Assumption 1.1. Let G{t) = g{s)ds. We assume 

(1) g{t) G C^{R+),g{t) > 0,g(0) < ^ and G{+oo) < . 

(2) There exists a convex function H{-) G C'^(M+), which is strictly 
increasing with H{0) = 0, such that 

g'{t) + H{g{t))<0,\^t>0. 

(3) Let y{t) he a solution of the following ODE 

y'{t) + H{y{t)) = 0,?/(0) = c/(0), 

and there exists oq G (0,1) such that G Li(]R+). 

(4) There exists 5 > 0 such that H{-) G (7^(0,5) and x‘^H"{x) — xH'{x) + 
H{x) > 0,Vx G [0,5]. 

(5) A satisfies ||u|| < Ao||A^/^u|| for all u & H. 

(6) 7 = a - 4^ > 0. 

Remark 1.2. Comparing to the result in [TB], here we do not require the 
convexity of g, however g{0) has to be suitably small. 

Theorem 1.3 (Existence of weak solution). Consider system ([I]). If 

{uo,ui,U2) G V{A^/^) X V{A^/^) X H, 

then, under the Assumption II.il this system has a unique weak solution u 
satisfying 

u G CHK+;E>(a 4) nC2(M+;17). 

Theorem 1.4 (General decay of energy). Under the Assumvtion \l.ll the 
energy of the weak solution decays to zero. Moreover, there exists positive 
constants C, 13, k such that 

E{t) < Cy{^t + k). 

1.2. Background. Moore-Gibson-Thompson(MGT) arises from modeling 
high amplitude sound waves. There have been quite a few works in this 
research field due to the wide range of applications such as the medical 
and industrial use of high intensity ultrasound in lithotripsy, thermother¬ 
apy, ultrasound cleaning, etc. The classical nonlinear acoustics models 
include Kuznetsov’s equation, the Westervelt equation and the Kokhlov- 
Zabolotskaya-Kuznetsov equation. A thorough study of the linearized mod¬ 
els is a good starting point for better understanding the well-posedness and 
asymptotic behaviors of the nonlinear models. Actually, the work m has 
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shown, even in the linear case, rich dynamics appear. In [19], Marchand et.al. 
2012 gave a detail analysis of this equation; using the abstract semigroup 
approach and a refined spectrum analysis they settle the well-posedness of 
m and identified an accumulation point of eigenvalues which essentially 
connects to the exponential decay of energy. Kaltenbacher et.al. m also 
studied the fully nonlinear version of MGT equation and gave the well- 
posedness and the exponential decay. 

In [ig, Kaltenbacher, Lasiecka and Marchand studied well-posedness and 
uniform decays of energy for the linearized MGT system, 

(3) TUttt + autt + c^Au -\- bAut = 0. 

2 

A critical parameter 7 = 0 ; — ^ was disclosed. It has been shown in m 
that when 7 > 0 , namely in the non-critical case, the problem is well-posed 
and its solution is exponentially stable; while 7 = 0 , the energy is conserved. 

It is well known that the wave equation conserves mechanical energy. 
More specifically, consider the wave equation 

utt — Art = 0. 

Multiply by ut, we have 

+ IIVup) = 0 ^ E(t) ^ i(||t.,|p + llVall") = E(0). 

Here we use E{t) to represent the mechanical energy, the summation of 
kinetic and potential energy. So the mechanical energy does not change 
when a wave evolves. On the other hand, we normally see mechanical energy 
dissipates in a physical system due to certain damping mechanism^ A 
different viewpoint is, in order to force the energy decay, we have to inject 
damping mechanisms into the system. In this aspect, there are different 
ways to implement the idea. One way is to add viscous damping, also called 
frictional damping, into the wave equation. Namely, taking into account the 
friction, we end up with equation in form of 

Utt - Au + Ut = 0. 

The energy estimate gives 

+ llVall") = -||u,||" < 0 ^ E{t) < E{0). 

Obviously, the mechanical energy does not increase for sure; in fact it decays, 
at least when the kinetic ||ui|P is not zero. Indeed, it can be shown, by 
standard Lyapunov function method, that the energy decays exponentially 
[ 8 |. Adding the structure damping Aut to the wave equation is another 
way to obtain exponential decay of the energy. These topics are of great 
importance in real applications and form an active research area [Tj. 

^Damping is the dissipation of energy, which transform the mechanical energy into 
another form, e.g. heat or light. 
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What interests us here is a different damping mechanism caused by vis¬ 
coelasticity, which forces the appearance of memory term in a system. Vis¬ 
coelasticity is the property of materials that exhibit both viscous and elastic 
characteristics when undergoing deformation. It usually appears in fluids 
with complex microstructure, such as polymers, suspensions, and granular 
materials. One encounters viscoelastic materials in biological science, mate¬ 
rials sciences as well as in many industrial processes, e.g., in the chemical, 
food, and oil industries. The phenomena and mathematical models for such 
materials are more varied and complex than those of pure elastic materials 
or those of pure Newtonian fluids, see [22]. 

In [18] , we investigate the case where memory effects are incorporated into 
the model. It is well known that memory generates stabilizing mechanism 
for second order wave equations. There is an abundant literature on the 
topic, see [2[ [S] IU [G] |5[ (9] jTOl jlU |20l [2T[ [23[ [2^, to name a few. The 
prototype model is viscoelastic wave equation, 


( 4 ) 


utt — Au -|- / g{t — s)Au{s)ds = 0, 

Jo 


where A is Laplacian defined on bounded domain with smooth boundary. 

The convolution term g{t — s)Au{s)ds represents memory: the integral 
itself suggests the nonlocality in time; the system at present moment “re¬ 
member” certain information in the past history. This is a simplified one 
dimensional wave equation with memory. Like all the physically meaningful 
partial differential equations, it is derived based on the balance laws of mo¬ 
mentum and/or conservation of mass, together with a constitutive equation 
relating stress to strain (elasticity) or stress to rate of strain (fluids). For 
more information on the derivation of this 1-D system, we refer to Renardy 
m- On a detailed description on damping mechanisms, see [T] and the 
references therein. 

In the study of wave equation, an interesting phenomena is that, when 
two different kind of damping terms appear simultaneously, despite the fact 
that both terms help the system dissipate the energy, they do not necessarily 
accelerate the decay process. For instance, while the frictional damping 
alone causes exponential decay; only polynomial decay can be reached when 
an additional memory term presents itself with a polynomial kernel]?]. 

Natural questions can be raised based on the study of wave-memory sys¬ 
tem: What kind of term creates damping effect in a system? How to select 
multipliers that can “detect” the damping? Does the damping term always 
help? Is it possible that adding a damping term “hurt” the initial energy, 
in any sense? While the answers to these questions in case of wave with 
frictional damping are relatively simple, they are not for the memory damp¬ 
ing we introduce below. Following the wave-memory system, this work is 
part of the effort in understanding this damping mechanism. To our best 
knowledge, this is one of the first few works that study the memory damping 
in the context of a third order in time system. 
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In the rest of this section, we recaii the resuits for MGT equation. In 
Section 2 we prove the main resuits stated before. 

1.3. Results for MGT. Consider the MGT equation ([3]), 

TUttt + (xutt + c^Au + bAut = 0, 

with initiai conditions 

U(0) = Uo,Ut(0) = Ul,Utt(0) = U2. 

Here M is a seif-adjoint positive operator on a Hiibert space H with a dense 
domain T){A) C H. The initiai data 

(uo,ui,U 2 ) £n = X H. 

We state the foiiowing resuits from m- 

Theorem 1.5 (Kaitenbacher et.ai. 2011, m)- Let T > 0, 6 > 0, a G 7^, 
then the system ^ generates a strongly continuous semigroups on Li. 

2 

Moreover, introduce the parameter 7 = 0 — we have 

Theorem 1.6 (Kaitenbacher et.ai. 2011, [I2])- Let r > 0, 6 > 0, a > 0, c > 
0. Let 

Ei{t) = h\\A^/‘^{ut + + T\\utt + 

E2{t) = a\\ut\\‘^ + 

and 

E{t) = Ei{t) + E2{t). 

1. Ifj > 0, then the semigroup generated in Theorem \1.5\ is exponentially 
stable on Li. And there exist a; > 0, C > 0 such that 

E{t) < Ce-^^E{<d),t > 0. 

2. If j = 0 , the energy Ei (t) remains constant. 

Remark 1.7. The MGT equation ^ can be rewritten as 

7^2 

TZtt + bAz + 'yzt = —r^t, 
b 

2 

where z{t) = ut + ^u. Since r > 0 ,6 > 0, it is easy to see that, in the new 

variable z, ^ becomes wave equation when 7 = 0 , hence no decay can be 

expected. This is exactly how MGT is connected to wave equation through 

2 

the parameter 7 = 0 — On the other hand, when 7 > 0, m becomes a 
damped wave equation which intends to force exponential decay. 

Remark 1.8. Theorems \ 1. 5\ and \1.6\ were firsted vroved in [12]. A different 
proof of Theorem \1.6i. which can serve a warm up for the current work, was 
given in |18j . 
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Remark 1.9. In this work, we work in the non-critieal regime (7 > 0). 
Both the “frictional damping” and “memory damping” are present in the 
system. It turns out that, while the frictional damping alone causes expo¬ 
nential decay, we are only able to get a slower decay if a weaker memory 
damping is added to the system. 

From now on, we study MGT with memory terms and prove the main 
results stated before. 


2. Proofs 

Remark 2.1 (On the existence, Theorem II. 3p . We omit the proof and fo¬ 
cus on the decay part. The concept of weak solutions and the corresponding 
functional setting are from |12] . The proof can be completed by the standard 
Galerkin method. The key is a global energy bound, which is a by product of 
the decay part proof below. 


Proof of the general decay, Theorem 11.41 

The proof is lengthy and quite involved. So we split the proof into three 
subsections. In the first subsection, we derive the energy inequality. All the 
later calculations are based on this estimate. However, we are not able to 
get decay rate out of this continuous version of energy estimate, for it is 
not in an form applicable to Gronwall inequality. To get what we want, we 
need two parts: 1. the discretized version of energy estimate developed in 
m, which essentially is a generalized version of Gronwall inequality; 2. the 
iteration technique developed in m, which can optimize the estimate in 
finite steps. In the second subsection, we discretize the estimate, based on 
which we are able to get an “initial” decay estimate. The decay rate here 
is not optimal. We improve the result in the third subsection through an 
iteration process. The iteration can be hnished in hnite steps and we are 
able to get the optimal decay rate. 

2.1. Energy inequality. In this subsection we carry out the energy esti¬ 
mate in several steps. The goal is to get a differential inequality in form 
of ^E{t) = —R{t), where E{t) is the natural energy functional and R{t) 
is a positive function called damper, it forces E{t) to decrease and creates 
energy damping. Generally, the expression of natural energy E{t) looks 
sloppy, hence standart energy funetional F{t) is defined and proved to be 
equivalent to E{t). Namely there are positive constants C'i,C'2 such that 
CiE{t) < F{t) < C 2 E{t). As a dressed up version of E{t), F{t) usually 
makes the calculations much neater. 


Recall the equation 

(5) TUttt + ciUtt + c^Au + bAut — / g{t — s)Au{s)ds = 0. 

Jo 




MOORE-GIBSON-THOMPSON EQUATION WITH MEMORY, PART II: GENERAL DECAY OF ENERGY 


Step I. Multiplying ([5]) by uu, we have 
d 


dt 


with 


Ei{t) = -2a||nft|p + 2c^\\A'-^‘^ut\\‘^ - 2g{t){Au,ut) 

+2 f g'{t - s){A{u{t) - u{s)),utit))ds, 

Jo 

Ei{t) = [r||nti|p + bWA^/^utW"^ + 2c^{Au,ut)] 


-2 


/ g{t - s){Au{s),ut{t))ds. 
Jo 


—E 2 (t) = -2b\\A^^‘^ut\\'^ + 2T\\utt\\‘^ + g' ° A^^’^u - g{t)\\A^/'^u\ 


Step II. Multiplying ([5]) by ut, we have 

A 

dt 

with 

E2{t) = + a\\ut\\‘^ + 2T{utt,Ut)] 

rt 


+goA^^'^u— f g{s)ds\\A^^'^u\ 

Jo 


Recall goA^^‘^u = g{t — s)\\A^/'^u{t) — A^^‘^u{s)\\'^ds. 

2 

Step III. Define natural energy functional. Since 7 = 0 ; — ^ > 0 is 
assumed, we can pick a k such that 

a 

— < k < -. 

b T 

Let E{t) = Ei{t) + kE 2 {t) be the natural energy., we have 


E{t) = h\\A^^‘^ut\\^+ 2(fi{Au,ut) + (?k\\jfi^‘^ 


u\ 


( 6 ) 


+ T\\utt + kut\\ +kT{ - k)\\ut\\ 

T 

+ kg o A^^'^u — kG{t)\\A^^‘^u\\‘^ 

+ 2 [ g{t — s){A^^'^{u{t) — u{s)),A^^‘^ut{t))ds 

Jo 

- 2G{t){Au{t),ut{t)). 


and 


^E(t) + R(t) = 0, with 
dt 

R{t) = 2{a - kT)\\utt\\'^ + 2{bk - c^)\\A^^‘^ut\\‘^ 
+2g{Au, Ut) + kg\\A^^‘^u \^ 

—2 f g'{t — s){A{u{t) — u{s)),ut{t)ds — kg' o 
Jo 
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Remark 2.2. The calculations appear to he lengthy. But two hints should 
make it less tedious: 1. since we are working on real Hilbert space, many cal¬ 
culations here are like algebras of completing square; 2. both E{t) and R{t) 
contain two parts, one part contains terms from the normal MGT equation, 
the other stands for memory terms. 

2 

Remark 2.3. Although the calculations here are for arbitrary k G 

we will need k to satisfy an additional condition stated later, see Lemma 

WMK 

Step IV. Define standard energy functional. Now denote the standard 
energy 

F{t) = ||u«||2 + + 9 o 

we claim that 
Lemma 2.4. E{t) ~ E{t). 

Remark 2.5. Here F{t) contains a memory term in form of g o A^^'^u, 
contrary to the energy functional defined in [?], which contains memory 
term —g' o Afl'^u instead. This is the consequence of us trying to eliminate 
the convexity restriction on g{t). 

Proof. We want to prove that one is bounded by a multiple of the other. 
E{t) = + 2c‘^{Au,ut) + c^k\\A^^‘^u\\‘^ 

+T\\utt + kutW"^ + kT{— - k)\\ut\\‘^ 

T 

+kg o A^^'^u — kG{t)\\A^^‘^u\\^ 

+2 f g{t — s){A^^‘^u{t) — A^^'^u{s),A^^‘^ut{t))ds 

Jo 

-2G{t){Au{t),ut[t)) 

= + kA^/M'^ + {b- ^l)\\A^Aut\\^ 

k k 

+T\\utt + kutW’^ + kT{- - k)\\ut\f 

T 

k 

+2 f g{t — s){A^^'^u{t) — A^^'^u{s),A^^'^ut{t))ds. 

Jo 

Remark 2.6. The summation of last three terms is non-negative, since 
2 

<2, f g{t — s)\/k\\A^^‘^u{t) — A^^'^u{s)\\^^\\A^^‘^ut{t)\\ds 

Jo \k 


f g{t — s){A^^'^u{t) — A^^‘^u{s),A^^'^ut{t))ds\ 

Jo 
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< [ g{t - s)[k\\A^^‘^u{t) - + Y\\A^/‘^ut{t)\f]ds 

Jo k 

= kg o AJ^'^u + with G{t) = f g{s)ds. 

k Jo 

Remark 2.7 (Energy match). While we can guarantee the non-negativeness 
of the summation, that is not suffieient for proving the equivalenee of E{t) 
and F{t). For that purpose, we expeet a positive eopy of g o AJI‘^u to come 
out. This can be done by energy match. The idea is to borrow a portion of 
from term (6 — ^)\\A^/‘^ut\\‘^ and add it to term ^^\\A^^‘^ut\\‘^. 
Then through completing square, we obtain a sub.stantial portional of g o 
A^/^u. 

Split the term (6 — ^)\\A^/‘^ut\\^ to get 

E(t) = + \{b - f )||A/V||" 

+T\\utt + kut\\^ + kT{— - A:)||ttt|p 

T 

+[)((.-j) + ^IIIAA.,IP + *9° 

+2 f g{t — s){A^^‘^u{t) — A^^'^u{s),A^^'^ut{t))ds. 

Jo 

Pick a do, so that 0 < ao < k and This can be 

done since G{t) is finite and we can choose ao sufficiently small. Then by 
Holder inequality we have 

2| [ g{t—s){A^^‘^u{t)—A^^‘^u{s),A^^'^ut{t))ds\ < {k—ao)goA^^'^u+ 

Jo \k — <To) 

So 

E(t) > + i(6 - d)||^i/2„,||2 

k 2 k 

+T\\utt + kutW^ + kT{- - k)\\ut\f 

T 

+[i((.-d) + M]||_4V2,„||2 + tj,o_4'/2,, 

-(fc - <r(,)goA^^‘^u - — 

[k - ao) 

> T:^||^i/2„, + ,,A'l\f + i(6 - d)||^i/2„,||2 

+T\\utt + feutip + kT{- - k)\\ut\\^ 

T 


+aog o A^^'^u. 
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Since (? — G{t) is strictly positive, we have E{t) > CF{t) by Lemma 12.81 
below. 

Lemma 2.8. Let Cq > 0. Then 

||/ + 5|P + Co|blP~||/|P + lblP. 

Proof. Note 

11/ + 5lP + C'ollfl'iP = ^ ll/ll^ + 2(/,5') + (1 + -^)lblP 

i + — ^ 

> Cidl/ll" + Hall"), Cl = mm{(l - ^}, 

The other direction \\f + g\\‘^ + C'oUfl'lP < ^2(11/11^ + Hfl'lP) is trivial. □ 
On the other hand, it is easy to see 

^ - gw 11^1/2^^ + + (b - 

k k 

+T\\utt + kutW"^ + kT{- - k)\\ut\\'^ 

T 

+9^\\AG\,\\‘^ + kgoA^/\ 

k 

+2 f g(t — s){A^^‘^u{t) — A^^‘^u{s),A^^‘^ut{t))ds 

Jo 

< C\\A^/\t + kA^/‘^u\\‘^ + C\\A^/\t\f 

TC'lliiii + 

+C\\A^/‘^ut\\^ + Cg o A^^‘^u 

< CWuttW^ + C\\ut\? + C\\A^l^ut\\^ 

+C||.4,^'^^rt||^ + Cg o A^^'^u 

< CF{t) 

because of the boundedness of g{t), G{t). The equivalence of E{t) and F{t) 
is established. □ 

Remark 2.9. Obviously F{t) is mueh simpler than E{t). Indeed, we have 
cleaned up the cross terms and the eoeffieients. Moreover, since ||t(;|| < 
Ao||.4.^/^t(;|| for all w ^ H, F{t) also absorbs terms ||ui|p and ||m|P. 

Step V. Analyze the damper. Recall 

±Eit) = -R(t) 

with 

R{f) = 2 ( q ! - kT)\\utt\\^ + 2{bk - c^)\\A^/‘^ut\\^ + 2g{Au,ut) 
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+kg\\A^^‘^u\\^ — 2 f g'(t — s){A{u{t) — u{s)),ut{t)ds — kg'o 

Jo 

We claim R{t) > 0 and hence E{t) is non-increasing. 

R{t) = 2{a - kT)\\utt\\'^ + 2{bk - c^)\\A^/'^ut\\'^ 

+ ^\\A^/^ut + kA^/\\\^ - ^\\A^/\t\\^ - kg' o A^/^u 

k k 

-2 [ g'{t - s){A{u{t) - u{s)),ut{t)ds. 

Jo 

For the last term, we can pick a number 5 G (0, k) such that 

-2 f g'{t - s){A{u{t) - u{s)),ut{t)ds 
Jo 

= 2 / — s)y/k — 5\/—g'{t — s)\/l/k — 6{A{u{t) — u{s)),ut{t)ds. 

Jo 

1 ft 

<—{k — 6)g^ o / g^{t — s)ds 

k — 0 Jq 

= -{k - 6)g' o A^/^u + 11^1/2^,||2 

k — 0 
So 

R{t) > 2(q: - kT)\\uu\\^ + ^^\\A^^‘^Ut + kA^/‘^u\\‘^ 

k — 0 

-5g' oA^/‘^u + {-j^^ - ^)g{t)\\A^/‘^ut\\^. 

Now, given the Assumption 11.11 by the Lemma 12.101 below, we can pick 
cr > 0 such that 2{bk — (? — becomes strictly positive. 

Lemma 2.10. //O < ^(0) < then there are a > 0 and k G (^, such 
that g{0) < {k — cr){bk — (?), or equivalently, 2{bk — c? — |^) > 0. 

Proof. Since g{d) < = y(6y — c^), all we need is to show 

{k — a){bk — c ) —(b - c ) as k ^ —,a ^ 0, 

T T T 

which is trivially true. □ 

Thus there exists a positive constant such that 

R{t) > Cs{\\utt\\^ + - g' o A^/'^u). 

Obviously, 

(8) - g'o A^/‘^u< ^R{t). 

This inequality will be used later, where we construct convex function and 
apply Jenson’s inequality. See lemmas [2.161 and 12.171 
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Moreover, integrate ([7|) on [t, T], we have E{T)+R{s)ds = E{t). Hence 
E{T) + Cs - g' o A^/^ds < E{t). 


In particular 


(9) 


||nn|P + - g' o A^^'^uds < CE{t). 


What we expect now is 




To achieve the goal, multiply ([I]) by u, we have 


dt 2 


u 


= a 


d ,T 


^ “ 'r{utt,u) - a{ut,u)] 


+G{t)\\A^^‘^u\\^ — f g{t — s){A^^'^{u{t) — u{s)),A^^‘^u{t))ds. 

Jo 

which gives 

{(? — G{t) — eG{t)) J 

rT ^ 

<ay \\ut\\‘^ds + [-\\ut\\‘^-T{utt,u) - a{ut,u)]\f 

By Assumption 11.11 1. G{t) < G(+oo) < cP, we can choose e so small that 
{P - G{t) - eG{t)) > C, > 0. 

Furthermore, 

/ WutW'^ds^Xl! \\A^/‘^ut\\'^ds <CE{t), 

Jt Jt 

\\utf\^ < \\A^/^ut{t)\\^ + \\A^/^ut{T)f < GE{t), 

\\{utt,u)\\ < llttttlp + ||n|p < GE{t), 

\\{ut,u)\\<CE{t), 

J g o A^^^uds < GE{t), 


hence 

( 10 ) 
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Now combine ([9]) and (jlOp . we have 

utt\\^ + — go A^^'^uds < CE{t). 

Remark 2.11. Consider the special case where H{s) = cs, hence g' < 
—eg, then we can get exponential decay without going further. Indeed, since 
—g' > eg, from the last inequality we have E{s)ds < C WunW"^ + 

— g' o A^^'^uds < CE{t). By a generalized Gronwall 
lemma, E{s)ds < CE{t) immediately gives the exponential decay. See 
Lemma 2.1 in m- 

Remark 2.12. In this work, we assume memory kernel g{t) has general 
decay rate, and we need techniques that is more sophisticate. To do the decay 
estimate, we adopt the idea in Lasiecka,etc [l6] and represent the decay by 
the solution of an ODE. For that reason we have to discretize the energy 
inequality and apply a series of lemmas. The technique has been applied in 

nil- 



step VI. 

The previous step implies 


So 


< CE{t). 

f E{s)ds <CiE{t) + C 2 f goA^^^^uds 

Jt Jt 

= CiE{T) + Cl / R{s)ds + C 2 [ go A^^^uds. 

Jt Jt 

Since E{t) is non-increasing function of t, we arrive at 
(11) (T - t - Ci)E{T) < Cl [ R{s)ds + C 2 [ goA^^uds. 

Jt Jt 

Remark 2.13. We emphasize the fact that the starting “moment” t is taken 
sufficiently large so that the energy estimate can be carried out. Other than 
that, t and T are arbitrary. The constants Ci , C 2 are independent of t and 
T. As we .shall see soon, the decay rate shall be discovered based on the 
relations between each pari of the following integrals 

f goA^Audg.!^ f [—gfoA^Auds^ [ R{s)ds. 

Jt Jt Jt 


2.2. Initial decay estimate. We introduce the Lemma 3.3 from m , which 
is the cornerstone of the later calculations. 
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Lemma 2.14 (Lasiecka & Tataru 1993, [M]). Let p be a positive, increasing 
function such that p{0) = 0. Since p is increasing, we can define an increas¬ 
ing function q,q{x) = x — {I + p)“^(x). Consider a sequence Sn of positive 
numbers which satisfies 

Sm+l + PiSm+l) < ■Sm- 

Then Sm < S{m) where S{t) is a solution of the differential equation 
±S{t) + q{S{t)) = 0,S{Q) = So. 

Moreover, if p{x) > 0 for x > 0 then limi_,.oo S{t) = 0. 

Remark 2.15. This lemma suggests that we are expecting an discretized 
energy inequality in form of E{Tn+i) + H{E(Tn+i)) < E(Tn), with an in¬ 
creasing non-negative function H. 

We put the calculations in several steps. 

Step 1: Discretize energy estimate 

In (fTT]l . pick t = Tn = nTo and T = T„_|_i = (n + l)To, where Tq is a given 
positive real number and n is any positive integer. Here Tq is chosen to be 
sufficiently large so that all the estimates before can be justified. We have 

Mn+l r^n+1 

(To - Ci)E{Tn+i) < Cl / R{s)ds + C 2 go A^/^ds. 

Here we pick Tq large enough so that Tq — Ci > 0 and we have 

fTn + l fTn+1 

(12) E{Tn+i)<Ci R{s)ds + C 2 goA^/^uds. 

J Tn '^Tn 

Note n is arbitrary positive integer and C'i,C'2 depend on Tq, but do not 
depend on n. 

Step 2: Construct convex function related to H 

For the oq £ (Oj 1) given in Assumption ll.il define function 

(13) Hi^aois) = aos^~^ Hafis) = aos^~^H{s^). 

Lemma 2.16. Denote ||.4.^/^?x(t) — A^/‘^u{s)\\ by f{t,s). Under the As¬ 
sumptions [T^ there exists an interval [0,5),0 < 5 < 5 on which we have 

1. Hi^aoi^) = 0 Hi^aois) is increasing and convex. 

2. Moreover, if Ca^ = supc(ao,t) = sup — s)f‘^{t, s)ds < 00 , 

t>0 t>0 

then there exist constants i? so that 

(14) Hi^ao['d{goA^u)]{t)<^^R{t),fort£[Tn,Tn+i], n = l,2..., 

^CT 

with 1? G (0,1) independent of n. 
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Proof. The proof was given in m- We present it here for readers’ con¬ 
venience. Since H E oo), it is easy to see that function = 

I _k j_ 

aos °‘o H{s°‘o)^0 < ao < 1 is well-defined on [0, oo) with = 0- 

1. Let k = -^ > 1. We have Lfi and 

= + [H'is’^) - s-^H{s^)] > 0,Vs > 0. 

K 

The last inequality results from the properties of H[s) 

(15) H{x) < xH'{x),x > 0. 


Indeed, this follows from geometric interpretation of convexity of H(x) and 
the fact that H{0) = 0 which then gives H'{x)x — H{x) > 0,Vx > 0. Thus 
we show that Hi^a^ix) is increasing on the positive half line. 

For the second derivative on (0, oo), we have 


+(1 - k)s-^H'{s'^)s'^-^ + ks^-^H''{s^) 

1 


(16) 


gk+l 

k-1 




[kx'^H"{x) — {k — 1)xH'{x) + {k — l)H{x)] 
{x^H"{x) - xH'ix) + H{x)) + -^H"{x), 


with X = s^. 

In view of Assumption ll.il we conclude the proof for part 1, namely, there 

I _L j_ 

exist an interval (0,(5),0 < <5 < 5 on which Hi^aois) = aos °‘oH{s°‘o) is 
increasing and convex. 

2 . 

I _L J_ 

Now we prove that Hi^aois) = aos °‘o H{s°‘o ) satisfies 


(17) 


o A 2 u){t)] < ^R{t) 
1 O'—1 ( ory& 


(or equivalently {g o A2u){t) < 00(7; '^(^)))’ under the assumption 

0 < CaQ = sup c(aoT) = sup / g^~°‘° {t — s)ds < 00 . 

t>o t>o Jo 

Recalling that {g o A 2 u){t) = g{t — s)ds and c(aoT) < 00, also 

1 -J_ J_ 

noting Hi^aois) = aos °‘oH{s°‘o) is convex for s small, pick tJ sufficiently 
small, by Jensen’s Inequality (see Proposition 13.11) we have 


H 


1,Q:o 


^ [ g{t - s)f{t,s)ds 

Jo 


= H 


1,0:0 


'dg'^°{t - s)g^ °‘°{s)f{t,s)ds 


= Hi, 


ao 


c{t^ O^q) 


[ 'dc{t,ao)g°‘°{t - s)g^ °‘°{t - s)f{t,s)di 

Jo 
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-f^i,ao ['&c{t,ao)g°'°{t - s)]c?^ °‘°{t - s)f{t,s)ds 


c{t 


c{t, ao) 

J ao['dc{t,ao)g°^°{t-s)Y~'^H{[dc{t,ao)g°‘°{t-s)]'^)g^~°‘°{t-s)f 


1 - 7 ^ rt 


act'd “0 


c{t,ao)°‘o -^0 


H{['&°‘oc{t,ao)°‘og{t - s)])f{t,s)ds 


Now we can make i? so small such that i?“o c{t, ao) “o < 1, thus we have 

c{t, ao)“° g{t — s)]) < '(^“o c{t, ao)°‘o H{g{t — s)), 
because H{0) = 0 and H(x) is convex. So we have 

'& [ g{t - s)f{t,s)ds <ao'd ! H[g{t - s)]f{t,s)ds 
^ Jo ^ Jo 


H 


1,00 


< 


ao'd 


ao^ j [-g'{t-s)]f{t,s)ds 


a 


The final inequality s from Q. This completes the proof of (jl7p . □ 

Lemma 2.17. Given the results of lemmas above, our energy functional 
E{t) satisfies 

( 18 ) ^(r„+l) + Hl,c,,{E{Tn+l)} < E{Tn) 

holds for all n > 0, where Hi^ao is defined by 




-1 

ao 


ao'd 


L^Tr 


aUO 


+ C'2X, Vx G R. 


Moreover, we can show that Hi^ag{s) is a convex, continuous increasing and 
zero at the origin function. Here Ci,C 2 ,'& depends only on ao,To, but not 
on n. 


Proof. From ()12p . we have 

E{Tn+l) < Cl [ 


'n+1 


{g o A' 2 u){ t)dt + C 2 


'n+l 


R{t)dt 


< 


Cl ^_i ,aoi? 

1 ? .L 


1 act'd f 

Hfil^i-^R{t))dt + C2 / R{t)dt 

J XU 


ao 




Here we define G by 


Lc-.r 




ao'd 

° / R{t)dt 

fTn+l 

/ R{t)dt 

'J Tn 


+ C 2 


in + 1 


R{t)dt 


Cl, 


r-l 


ao'd 
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It is easy to see that G is increasing and concave near 0 with G'(O) = 
0. Denote G~^ by which is increasing, convex and through origin. 

Furthermore, 


E{Tn+i) < G 


Hl,ao[E{Tn+l)] < 


ira+l 


R{t)dt 


cT, 


n+1 


R{t)dt 


Hi,^,[E{Tn+l)] < E{Tn) - E{Tn+l). 
E{Tn+l) + Hl,o.o[E{Tn+l)] < E{Tn). 


□ 


Step 3: Initial estimate 


Define H„,^ =/—(/ + Hi^ao) By Lemma 12.141 and Lemma 12.171 we 
immediately have 

Lemma 2.18. Given the lemma \2.n\ and condition we have the fol¬ 
lowing decay rates for the energy function 

E{t) < s(t),Vt > T 

with 


(19) St+ ^Q,jj(s) — 0,s(0) — i?(0). 

Moreover, we have 


Lemma 2.19. All three convex functions, (x). Hi^aoix) and Hi^aoix), 
have the same end behavior at origin. 


Proof. Note the following 
H^,=I-iI + = (^ + ^l,ao) ° 

O (1 + Hl^ao) 

Since Hi^ao S ^^^[OjOo), we have Hi^ao{x) = 0{x) at the origin. We use 
A ^ B to represent that A and B have the same end behaviors. Then 

(/ + Hl^ao) ~ ^ ^ Eao O (/ + Ilceo) ~ ® ^ ^ ~ E\ aQi X ^ 0. 

For a detailed discussion, see [151 Corollary 1. p. 1770]. 

By similar arguments, from the relation between Hi^ag(x) and Eli^aoix) 


H 




-1 

«o 


at} 


G„T< 


0 J 


+ G 2 X, 


and the fact that they are both convex at zero, we can see ~ 

Ei^aoix) for X > 0 in the neighborhood of 0. This complete Lemma [2.19[ □ 


Corollary 2.20. There zs do > 0, /3 > 0 such that when 0 < x < do, 

Kaoi^) > PHi^aoix). 
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2.3. Improve the decay rate. From Lemmas 12.181 and 12.191 we know how 
the natural energy decays: it is bounded by a decay function s{t), which is 
the solution of an ODE. The ODE portraits the decay property of its solution 
through the convex function Hi^ao{x), which has the same end behavior as 
Hi^aoix) near origin. So the decay rate is characterized by Hi^aoix)- While 
this information is quite useful, it is not optimal. We expect the decay rate 
to be characterized by H{x), which has a better decay than Hi^aoix) unless 
ao = 0. 

Our next step is to improve the decay rates through an iteration process. 
Eor that reason, we need a comparison lemma. 

Lemma 2.21 (Comparison lemma). Given y{t) satisfying 

(20) y'{t) + H{y{t)) = 0, y(0) = yo > 0, 
and s{t) satisfying 

(21) s'{t) + H^a{s{t)) = 0,s(0) = So > 0,Q!o G (0,1]. 

Both function y{t) and s{t) are positive, decreasing. Then we have 

(22) s{t) <y°°{l5t +k), 

for some constant f3,K, meaning the decay rate of s(t) is identical to y°‘°{t) 
up to an affine transformation in the coordinates system, which does not 
change the long time behavior. 

Proof. Let TL{y) = > 0. It is easy to see TL is positive and 

decreasing function on [0,oo) with dTL{y) = —Thus 

yt + H{y) = 0 ^ —$r\ = dt 

d'H{y) = dt 
=> niy) -TL{yo) = t, 

=> y{t) = n~^{niyo) +t). 


St + Kaoi^) — 0 
^ St ^ /dHi^aois) 

=> {s^)t < -/dH{s'^) 



-i/(s“o) 


1 J- 

^^(s“0) -^(s“0) > pt 

^s^{t)<n-\n{s^) + pt) = yipt + K) 
=> s{t) < y"°(/3t + K), 


Similarly, 
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with K = ?^(sg° ) — ^(yo)- C 

Apply Lemma [20] with oq = 1, we have 

Corollary 2.22. In view of the assumptions on g{t) and y{t), we have 
9 {t) < y{t). 

The previous lemmas show that the energy is bounded by a function s{t), 
which in turn is bounded by function y°‘°{t) as t ^ oo. It might contain 
a delay, but it makes no difference since we are considering the asymptotic 
behavior. 

To conclude the proof of Theorem 11.41 we adopt an iteration process. 

Lemma 2.23 (Iteration for optimality). In finite steps, we are able to get 
the following deeay rates for the energy function: There exists positive con¬ 
stants C, fi, K such that 

E{t) < Cy{j3t + R). 

Proof. By Assumption 11.11 and Corollary 12.221 we have 

c{ao,t) = f g^~°‘°{t - s)f{t,s)ds 
Jo 

<2E{0) [ g^-^°{t- s)ds = 2E{0) [ g^-^°{s)ds 
Jo Jo 

poo 

<2E{0) y^-^°{t)dt <oo. 

Jo 

This is the critical estimate to initialize the second part of Lemma 12.161 
Then by Lemmas 12.17112.211 we get 


(23) E{t)<s{t) = Ciy^°il3it + Ki). 

So we have the first decay estimate. To go to the second iteration and apply 
Lemma 12.161 with L7i^2oo’ need 

c{2ao,t) = f g^~‘^^°{t — s)f'^{t,s)ds < oo. 

Jo 

This is true since, by (f23]l 

c{2ao,t) = f g^~^°‘°{t-s)f{t,s)ds 
Jo 

<2 f g^~‘^°‘°{t — s)E{s)ds 
Jo 

< 2Co [ g^-^^°{t - s)y^°{/3is + Ki)ds < oo. 

Jo _ 

The last inequality is from Proposition 13.31 in the Appendix. So we can 
apply Lemmas l2.16N2.23] aerain to get E{t) < € 29 "^^°{ht + ^2)- 
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We can continue this iteration with ak = kao, until k fist reaches a 
number m such that am < 1 < cim+i- This means 

E{t)<Cmy^^°{Pmt + Km), 

hence ^ ^ 

sup / E{s)ds < Cm sup / + Km)ds < oo, 

t >0 Jo t>o Jo 

which is guaranteed by Proposition 13.31 It follows that 


c(l,t)=sup f s)ds < sup f 2E{s)d. 
t>o Jo t>o Jo 


s < oo, 


in the next iteration when we apply Lemma l2.16[ we simply pick Hi^i{s) = 
H{s) to start the final iteration. And we end up with E(t) < Cy{/3t + k). 

□ 


3. Appendix 

Proposition 3.1 (Jensen’s inequality ). Let LI he a measurable subset of 
the real line. Let f{x) be a non-negative function on with k = f{x)ds 
finite. If g{x) is a real-valued measurable function on and function ip is 
convex on the range of g. Then we have 

X 5(a:)/(x)dx) < ^ T{ 9 {x))f{x)dx. 

Remark 3.2. Search Jensen’s inequality online for the proof. 

Proposition 3.3 (a-Sequence). Let oq £ (0) 1) hit) £ C[0, oo) be a 
positive function decreasing to zero. Moreover 

y^-^°{t)dt = L < oo. 

Let m he a positive integer such that mao < 1 o.nd (m + l)ao > 1- Tet /3, k 
he positive numbers with f3 > 1. A finite sequence of functions, in form of 
definite integrals on [0,t], are generated in the following way: 

(24) 4(t) = [ y^-^^°{t-s)y^'^-^^^°{f3s + K)ds,k = l,2,...,m, 

Jo 

(25) Im+i{t) = f y'’’^°‘°{lJs + K)ds. 

Jo 

Then each Ikit), k = 1, ...,m + 1, is bounded, uniformly in t, namely, 

sup Ik{t) < oo. 
t>o 

Proof. We first assume y{0) < 1. 

For k = we have 

h = f {t - iPs + K)ds 

Jo 
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ri/2 


yl-ka,^^_^^y{k-l)c.^^l^s + K)ds+ f + 

Jt/2 


ftl2 ft 

JO Jt/2 

rV2 ft 


ft/^ ft 

/ 2/i-“o(s)ds+ / s)(is 

JO Jt/2 


= 2 / y “°(s)ds < 2L < oo. 

JO 

Above we use the fact that y(t) is a decreasing function, hence y{j3s + k) < 
vis). 

For Im+i, we have 

Im+i = [ y'^°‘°{/3s + K)ds < f y^~°‘°{/3s + K)ds < L < oo, 

Jo Jo 

since (m + 1)q;o > 1 ^ mao > 1 — ao and y{t) < y{0) < 1. 

For y with y{0) > 1, we can always find a time to so that y{t) < 1 on 
[to, oo), since y decreases to zero. Whether or not the integrals are finite only 
depends upon the asymptotic behavior of y at infinity. So y being larger 
than 1 on a finite interval [0, to] does not bear influence on our result. □ 


The research of L Lasiecka and X. Wang has been partially supported by 
NSF Grant; DMS 0104305. 
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